Polydisperse chains placed on a one-dimensional lattice 
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Using a transfer matrix technique, we calculate the entropy of polydisperse chains placed on 
an one-dimensional lattice, as a function of the density of internal and endpoint monomers. The 
polydispersivity is determined considering different activities for the two types of monomers, as is 
usual for equilibrium polymerization and living polymers. If the mean number of monomers in the 
chains is fixed, the entropy as a function of the density of monomers displays a maximum and is 
limited from below by the entropy of monodisperse chains. The increase of entropy due to the 
polydispersivity is a linear function of the density of monomers. We also calculate the distribution 
of chain sizes, obtaining an exponential relation. 

PACS numbers: 65.50.+m,05.20.-y 
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I. INTRODUCTION 

In the microcanonical ensemble, the entropy of a sys- 
tem may be calculated by evaluating the number of mi- 
croscopic configurations as a function of the internal en- 
ergy and other extensive variables, such as the volume 
and the number of particles. The entropy, expressed in 
terms of these variables, is a fundamental equation of the 
system and thus all its thermodynamic properties may 
be obtained from this equation. A particular problem of 
this kind is the study of the thermodynamic properties 
of flexible chains placed on regular lattices. This prob- 
lem has a rather long history, which reaches back to the 
thirties, where the thermodynamics behavior of diatomic 
molecules (dimers) adsorbed on surfaces was considered 
by Fowler and Rushbrooke p] . In the simplest version of 
the model, only excluded volume interactions are taken 
into account, and thus each site of the lattice may be 
either empty or occupied by a single monomer. In this 
case, the internal energy of the system will be constant 
and the temperature is not defined. The dimer prob- 
lem may be generalized by considering chains with more 
than two monomers. More precisely, we may define the 
entropy of a system of chains with M monomers each 
(M-mers) placed on a lattice with V sites as 



(M = 2) and fully covered (p = 1) two dimensional 
lattices, the entropy was calculated exactly 0, in other 
cases estimates were found through series expansions Q , 
closed form approximations 4] and transfer matrix cal- 
culations with finite size extrapolations p|. 

Here we will study the one-dimensional version of the 
problem. After discussing briefly the monodisperse case, 
where all chains are composed by precisely M monomers, 
we include the possibility of polydispersivity, allowing for 
a distribution of different numbers of monomers in chains. 
In particular, we consider an ensemble of polydisperse 
chains determined by different activities for endpoint- 
and internal monomers of the chains, such as is usual 
in models of equilibrium polymerization and living poly- 
mers |y, |3 • Perhaps the most studied experimental real- 
ization of equilibrium polymerization is the phase tran- 
sition in liquid sulphur, and a good review of this and 
other experimental systems is due to Greer |8|- 

In the following section we will define the problem in 
more detail and show how it may be solved in the micro- 
canonical ensemble. Section ITTT1 presents the calculation 
in the grand-canonical ensemble, using a transfer matrix. 
Besides obtaining the entropy, we calculate also the dis- 
tribution of chain sizes. Final discussions may be found 
in the conclusion. 
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where T(N p ,M, V) is the number of ways to place N p 
chains with M monomers in each on the lattice and the 
thermodynamic limit is taken with fixed density of occu- 
pied sites p — N p M/V. In the particular case of dimers 
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II. DEFINITION AND MICROCANONICAL 
SOLUTION OF THE MODEL 



Let us start with the formulation and solution of the 
monodisperse case. In this case we must calculate the 
number of ways to place N p chains with M monomers 
each on a one-dimensional lattice with L sites. Since in 
this case there are L — MN p empty sites, it is easy to 
conclude that this number is equal to 
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Taking the thermodynamic limit [L — > oo) with the frac- 
tion of occupied sites p = MN p /L kept fixed and using 
Stirling's asymptotic form, we get 
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(3) 



In Fig. ^some curves of the entropy as a function of the 
density p are shown. The entropy vanishes both at p = 
and p = 1 and has a maximum at a density p m which 
grows monotonically with the molecular weight M, as 
may be also seen in the figure. In the polymer limit M 
oo we find p m = 1 , but the entropy vanishes identically. 




FIG. 1: a) Entropy as a function of the fraction of occupied 
sites for monodisperse chains with molecular weight M be- 
tween 1 and 10. b)Density of maximum entropy as a function 
ofl-l/Af. 

The entropy per site as a function of the density is a 
fundamental equation of the system, and all thermody- 
namic properties of this gas of M-mers may be obtained 
from it. Being athermal, the temperature is not defined 
for this gas, but the intensive entropic variables (f/T 
and p/T) may be obtained taking partial derivatives of 
sm(p)- As an illustration, we will obtain the mechani- 
cal equation of state explicitly. Considering the lattice 
constant to be equal to a, the physical entropy will be 
S = ksaLsMip), where fcs is the Boltzmann constant. 



Since the system is one-dimensional, the pressure will be 
replaced by the force / applied to the chain, therefore 
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The result is 
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It is interesting to expand the right side of this equation 
in powers of p, and the first terms are 
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The first term corresponds to the ideal gas result, being 
equal to the number of M-mers divided by L. The higher 
order terms are all positive. In Fig. [21 the equations of 
state for dimers and tetramers are plotted. The linear 
behavior at low densities and the divergence as p — > 1 
are apparent. 




FIG. 2: Mechanical equation of state for dimers (full line) and 
tetramers (dashed line). 

In equilibrium polymerization or living polymers, the 
chains grow through a process in two steps: first a chain 
is started and then it grows incorporating an additional 
monomer each time. In the case of liquid sulphur, the 
first step corresponds to the opening of a small ring if 
atoms (mainly S s ), followed by the opening and addition 
of other rings to the linear chain. This process may be 
parametrized by two activities or chemical equilibrium 
constants. The first one, z' = exp(// /kgT), is associ- 
ated to the initial process and will be considered the sta- 
tistical weight of chain endpoints. The second activity, 
z — exp(p/ksT), relates to the addition of a monomer 
to an existing chain, and will be the statistical weight 
of internal monomers of the chains, p and p' are the 
chemical potentials of internal and endpoint monomers, 
respectively 6] . Usually the propagating reaction which 
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implies growing of the chains is faster than the opposite 
depropagating reaction, and this allows for the possibil- 
ity of synthesizing chains with a very narrow distribution 
of molecular weights by artificially terminating the poly- 
merization process, through the addition of an appropri- 
ate monomer [8j. If the equilibrium situation is reached, 
however, a polydisperse set of chains will be attained, 
and this is described by the simple model we consider 
with two activities. 

In the particular case of sulphur, the endpoint 
monomer activity z' has a very small value, of the or- 
der of 10 -12 , and therefore very long chains are formed 
in this polymerization process. The equilibrium polymer- 
ization model may be mapped onto the n-vector model 
of magnetism, and the activity z' will be proportional 
to the magnetic field of the n-vector model in the limit 
n — > 0. Thus, in this case the polymerization is a contin- 
uous phase transition |fg. For finite values of z', where 
finite chains are dominant, no phase transition occurs, as 
we expect for a lattice gas of chains with excluded volume 
interactions only. 

Now let us return to the calculation of the entropy 
in the one-dimensional model of equilibrium polymeriza- 
tion. The densities of endpoint monomers x and internal 
monomers y are conjugated to the activities z' and z, re- 
spectively. All configurations with the same numbers of 
chains N p = Lx/2 and internal monomers JVj = Ly have 
the same statistical weight. We may then calculate the 
number of these configurations in a combinatorial way. 
Besides the multiplicity in Eq. |J2J, an additional factor 
corresponding to the number of ways to place the internal 
monomers in the N p chains is now present. The result is 
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In the thermodynamic limit, we find the following ex- 
pression for the entropy per site 

s{x,y) = (1 - x/2 - y) ln(l - x/2 - y) - 

(1 — x — y) ln(l — x — y) — x ln(x/2) — 
ylny+(x/2 + y)ln(x/2 + y). (8) 

The mean number of monomers per chain will be M = 
2 + 2y/x, and if this number is fixed, we may express 
the entropy in the polydisperse case as a function of the 
density of monomers p = x + y. The result is 
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As expected, the entropy for the polydisperse case is 
never smaller than the one for the monodisperse case 



when M = M . The difference between these entropies is 

*M~*M = ^[(M-l)hi(M-l) 

-(M-2)ln(M-2)]. (10) 

The origin of the additional entropy in the polydisperse 
case is the second factor in expression 

Another point in this problem which may be addressed 
is the distribution of sizes of the chains. Again it is possi- 
ble to answer this in a combinatorial way, if we consider 
the probability that a particular chain has k = M — 2 
internal monomers. This probability corresponds to the 
ratio between the number of internal monomer configu- 
rations when one of the N p chains has exactly k internal 
monomers and the total number of configurations, which 
corresponds to the second factor in Eq. Q . In the ther- 
modynamic limit, this leads to the result 
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x + 2y \_x + 2y 
with k = 0, 1, 2, . . .. This result may be rewritten as 
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(11) 



where M = 2, 3,4, . . . and Tm is the fraction of chains 
with M monomers, which is a function of the mean num- 
ber of monomers in chains M. It is worth to notice that 
the Flory-Huggins approximation to a lattice model for 
equilibrium polymerization also results in an exponential 
distribution of the chain sizes [7| . Both results are in fact 
equivalent, if one takes into account that Dudovicz et al. 
considered isolated sites to be chains with one monomer, 
while here we called them empty sites. 



III. TRANSFER MATRIX SOLUTION IN THE 
GRAND-CANONICAL ENSEMBLE 

The grand-canonical solution of the monodisperse case 
may be found in reference p|, so the polydisperse case 
will be discussed here. The equilibrium polymerization 
model in one dimension was discussed already in refer- 
ence 0, with emphasis on the limit of infinite chains, 
where a phase transition occurs. Here we will use a sim- 
ilar transfer matrix technique to obtain the entropy and 
the distribution of chain sizes. To define a transfer ma- 
trix, we may use lattice gas variables rji = 0,1 associated 
to the bonds of the lattice, such that rji = 1 if the bond 
i is inside a chain and rji — otherwise. This definition 
is illustrated in Fig. [31 where a section of the lattice is 
depicted. 

The transfer matrix for this case is given by 



T = 
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FIG. 3: Section of the lattice with a particular configuration of 
chains. The values of the bond variable fi are indicated above 
the line and below the line the activity of each monomer is 
given. 



FIG. 4: Section of the lattice with a particular configuration 
of chains. The bond variables rj and monomer activities are 
shown, defined in a way to allow for the calculation of the size 
distribution of chains. 



The largest eigenvalue of this matrix is 



Ai 



l + z + v /(l-z) 2 + 4(z') 2 



(13) 



This eigenvalue may become degenerated when z 1 van- 
ishes, originating the polymerization transition studied 
in reference 9]. For our present interest, we will consider 
finite values of z' ', so that no degeneracy occurs. In the 
thermodynamic limit, the grand-canonical potential will 
be given by 



$ = -k B TVhxXi, 



(14) 



and since the entropy S is the partial derivative of <f> 
with respect to the temperature, and remembering that 
z = exp(/i/fcsT) and z' = exp(///fceT), we obtain the 
following expression for the adimensional entropy per site 



S z dX 1 z' d\i 

j—^ - lnAi - —^-q^ - ~ i ~q^J- ( 15 ) 



Now we may recognize that 
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z dXi 

y = t — s- 

Ai oz 

is the density of internal monomers and 
X ~ Ar dz> 

is the density of endpoint monomers in chains. Some 
algebra will lead us to expression [S] for the entropy in the 
polydisperse case. In the limit y — > we have x = p and 
the entropy of dimers sm=2(/°) is recovered. 

To obtain the distribution of internal monomers among 
the chains, we may define an activities for the internal 
monomers which are dependent on the localization of 
the monomer in the chain. Thus, the leftmost internal 
monomer will have an activity z% , his neighbor to the 
right has an activity Z2 and so on. The bonds of the lat- 
tice will also be numbered accordingly. A bond which is 
not inside a chain is represented by 77 = 0, as before, but 
the j'th internal bond, counted from the left to the right, 
is associated to r\ = i. These definitions are illustrated 
in Fig. |2 

The transfer matrix for the model now has an unlim- 
ited size, but a very simple structure 
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/ 1 z' 
z' zi 
z' z 2 
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(18) 



The secular equation may be found developing the deter- 
minant \T — XI I recursively by the first line. It is more 
convenient to express the equation in terms of to = 1/A, 
and after some algebra, we get 



-i+ w +(^)'[i+£nv) =0. (19) 

i=i j=i 

Now the density of monomers in the fc'th position in 
chains will be 
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X dz k lo dz k 



(20) 



Deriving the secular equation with respect to z k and then 
making all z k equal to z, we get 
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(z')M 1 - uz) 



{I - LJZ) 2 + {z') 2 L0{2 - LOZ) 



(coz) k . (21) 



This result may be written as a function of the monomer 
densities x and y 



Pk = 



x ( 2y 



2 \x + 2y 



(22) 



From this expression, the fraction r k of chains with k 
internal monomers expressed in Eq. Ill 1H may easily be 
found. 



IV. CONCLUSION 

Although it was rather simple to find the entropy and 
distribution of chain sizes in the equilibrium polymer- 
ization model using combinatorial arguments in the mi- 
crocanonical ensemble, it is interesting to perform these 
calculations also in the grand-canonical ensemble, since 
the combinatorial calculations are difficult to generalize 
to more complex situations, such as the problem defined 
on strips of finite widths (ladders). For the monodis- 
perse case, such calculations have lead to rather pre- 
cise estimates in the two-dimensional limit 0- We are 
presently performing similar calculations for directed and 
self-avoiding chains on strips, which hopefully will allow 
us to estimate the properties of these models in the two- 
dimensional limit. 
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